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ABSTRACT 



A probabilistic event-step computer simulation of a 
repairable item inventory system with exponential inner- 
arrival times between demands is constructed. The model 
allows for a wide variety of repair times and lead times 
to be considered. The major parameters are investigated 
for their sensitivity to system changes in order to help 
the user evaluate given or proposed inventory policies 
and parameters. A complete description of the model is 
presented and various measures of supply performance such 
as the expected number of unit years of backorders per 
unit period and the fill rate are provided. All flow 
charts and the GPSS/360 program listing are included. 
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I. INTRODUCTION 



Traditional inventory analysis is ultimately concerned 
with two decisions: (1) when to order and (2) how much to 

order. Most of the work done in the past addresses itself 
to these -questions for non-repairable items (consumables). 

A thorough treatment of such inventory models for both the 
backorders and lost sales case for numerous resupply and 
demand distributions now exists. 

However, at least in military applications, repairable 
items have now been found to be of great importance both 
materially and financially. Repairables often account for 
a significant portion of the items controlled by an activity 
and even more significantly when their portion of the total 
dollar investment is considered. For example, focus on one 
of the inventory control points of the U.S. Navy dealing 
with repairables such as the Aviation Supply Office (ASO) 
in Philadelphia, Pennsylvania. At fiscal year's end 1971 
repairables represented 18.7% of the different types of line 
items controlled by this activity for a total of 66,000 
separate items. These same items reflect an investment of 
nearly 3.1 billion dollars, almost 60% of ASO's total dollar 
investment. Further, ASO spends about 16 million dollars 
per month for repairs on these items and another 28 million 
dollars monthly for procurement of such items. Among ASO 
managed repairable items, 365,000 units were reworked in 
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fiscal year 1971. In view of the above figures it is easy 
to see why a need exists to help the manager improve his 
control over repairable items. It appears that an effective 
control over repairs would free for other needs many ASO 
dollars now spent on procurement. 

Until fairly recently the control of repairable items had 
not received much attention. Even in the case of consumables 
there are major difficulties in applying inventory theory to 
the problem because of the intricate relationships of the 
involved item parameters, such as demand and lead times. 

Also the number of different types of line items complicates 
the task. The introduction of repairable items makes the 
problem even more difficult since repairable items often 
require special storage, handling, packaging and shipping 
procedures due to their complexity. In addition, the 
introduction of the extra source of supply, the repair 
department, causes further complications. On the other hand, 
repairable items tend to minimize the number of different 
type of line items needed to be carried by an activity; the 
parameters are usually more easily determined than those of 
consumables; and the procurement lead times are more reliable. 
In addition, upper bounds for the total demand can often be 
determined, for if a given repairable item is a subcomponent 
of a larger unit, the number of larger units in use is 
generally easy to determine. 

Few people have addressed the problem of repairables. 
However, Schrady [Ref. 8] and others at the U.S. Naval 
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Postgraduate School, as well as Feeney and Sherbrooke [Ref. 6], 
have looked at the problem for the deterministic case. 

Feeney 1 s and Sherbrooke's model included the capability of 
handling compound Poisson demands but without degradation of 
supply. Wilson and Richards [Ref. 10] extended this model 
to allow for this possibility. Due to problems of convexity 
and others, no one has come up with a quick and easy way to 
evaluate and compare alternative repairable item inventory 
systems on a larger scale. Richards [Ref. 7] addresses this 
problem and this paper supplements his work by providing a 
detailed simulation model along with a programming guide and 
flow charts. This model enables the user faced with such 
problems to easily come up with numerical values for various 
measures of performance such as the average number of back- 
orders per unit period or the average number of fills per 
unit period. These measures of performance can then be used 
by the manager to evaluate and compare alternative systems. 
Further, the model can also be used as a teaching aide to 
help students and/or supply personnel get a clearer under- 
standing of the relationships between the various parameters 
involved in various inventory systems. A more detailed 
description of the model and how to use it are given in 
sections II, V, and VI of this paper. 
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II. THE MODEL 



The model presented in this paper is a general event-step 
computer simulation of the interactions between the major 
parameters for Q R, continuous review inventory systems. 

Here, Q is the order quantity and R is the reorder point for 
a given system. The purpose of this general model is to test 
and investigate the sensitivity of the parameters involved 
to help the decision maker evaluate, modify, and/or set up 
appropriate inventory systems. It is also used to demonstrate 
the analytical results obtained by Richards [Ref. 7]. The 
inventory systems simulated are single item, single echelon, 
and repairable or non-repairable systems where failures are 
assumed to be Poisson distributed. The model is a probabil- 
istic simulation in that the success or failure of any 
probabilistic event is determined in the model by comparing 
the numerical value assigned to the probability of success 
or failure to a program generated random number. In this 
model the numbers so used are obtained from many different 
distributions, using the one which best fits the particular 
parameter being considered in each case. Further, it is an 
event-step model in that the master program is controlled by 
a subprogram called an event list which is repeatedly updated 
as the master program switches from one event subroutine to 
another. This eliminates the necessity to look at the status 
of the entire system after each time unit (one hour) to 
decide whether or not an event has occurred. 
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Figure 1 shows the general type of inventory systems 
addressed in this paper . In order to help the reader under- 
stand the relationships pictured there, the following 
definitions for the on-hand inventory, repair inventory, 
and others are given. Note that because of the addition of 
the repair facility, the commonly used definitions of the 
various states for the consumable item inventory system 
must be modified. 

Definition 2.1 : The on-hand inventory (D-B) includes 

all of the items ready for issue (D) less any backorders (B) . 
Note that the on-hand inventory will be positive if there are 
no items on backorder. 

Definition 2.2 : The repair inventory (C) includes all 

items at the repair facility, including both those undergoing 
repair and those waiting to begin repair. 

Definition 2.3 : The net inventory (D-B+C) includes the 

on-hand inventory and the repair inventory. 

Definition 2.4 : The inventory position (D-B+C+A) is 

the sum of the net inventory plus all items which have been 
placed on order but have not yet been received. 

Consider the inventory system as one which has control 
of a large number of items of a single type where the 
number of demands (failures) for any time period is a value 
of a random variable having a Poisson distribution with 
mean rate X* The items which fail are assumed to be 
repairable with probability P and with probability (1-P) 
the items are scrapped. These failures will be referred to 
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as type one and type two, respectively. The manager 
controlling the system attempts to minimize the time 
required to fill a demand by maintaining a stock of items 
ready for issue. Systematically, the manager supplements 
his stock by procurement of new items and/or receipt of 
repaired items. Hence the manager is faced with the problems 
of when and how much to reorder. Standard approaches try 
to answer these questions by concentrating their efforts on 
minimizing the total expected cost. However, a difficulty 
in this approach is the determination of proper weights to 
assign to the various cost factors involved. 

A second approach concentrates on a class of simple 
procurement policies based on one or two parameters. The 
goal of this approach is to describe the long-term effect 
of an arbitrary member of the class of simple procurement 
policies being considered and then determine the choices 
for the corresponding parameters which are best with respect 
to a given measure of performance. This computer simulation 
allows the user to do this quickly and easily. 

After a class of simple policies has been specified, 
the study of the inventory model becomes a study of the 
associated stochastic processes and the long run properties 
of the processes. Emphasis is given to determination of 
stationary distributions which describe the probabilities 
of being in the various states after the system has been in 
operation for a long period of time. This computer program 
greatly simplifies the procedure whereby actual numbers are 
produced for these various entities. 
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The program is capable of representing the repair 
department in one of two ways. First, as an infinite- 
server queue with general independent repair times which 
can be constant or realizations of a uniform, normal, or 
empirical distribution. Under this policy the repair 
department acts as a M/G/«© queueing system which for the 
random repair times case will allow crossing of repair times. 
The crossing of orders is also allowed in this case. 

The second method considers the repair department as a 
single server queueing system. That is, only one batch of 
items up to a maximum quanitity (K) can be repaired at one 
time. Under this policy the user also has the option to 
enforce a constraint on the minimum batch size (L) as* well. 
This latter policy is provided because frequently in the 
case of repair, large set-up costs, tool-up costs, or 
manpower constraints make the second repair policy more 
economical . 
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III. EXPERIMENTS 



A. OBJECTIVES AND APPROACH 

1 . Obj ectives 

The objectives of this paper are threefold. (1) To 
show the need for "and justification of such a simulation. 

(2) To demonstrate the analytical results obtained by 
Richards [Ref. 7] . (3) To provide a means through which a 

decision maker can evaluate given or hypothetical inventory 
systems . 

The first objective has been touched on already and 
will be complemented by the rest of this paper. Consistent 
with the second and third objectives, a number of experiments 
were performed which involved generating operating statistics 
in order to quantitatively compare the simulated results to 
the theoretical and display the sensitivity of various 
parameters for particular systems. 

2 . Approach 

Each experiment consisted of a series of simulation 
runs. Each run simulated up to 14 unit periods, where a 
unit period represented ten years of actual operation of a 
given inventory system. For each run, operating statistics 
were collected for all measures of performance listed in 
Table I. For illustration purposes two measures of per- 
formance were chosen to judge the performance of the hypo- 
thetical systems examined in this paper. Those chosen were 
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TABLE I. SUPPLY PERFORMANCE MEASURES 



CF) = XH100: The average number of fills in a 

unit period. 

(FILRT) = XH99 : The fill rate multiplied by 100. 

This is the probability of being able to satisfy a 
demand without a time delay. 

(B) = XH98 : The average number of backorders 

incurred in a unit period. 

(Bj) = XH97 : The expected number of unit years 

of backorders in a unit period multiplied by 100. 

(BARB) = XH96 : The average number of items on 

backorder at an arbitrary point in time multiplied 
by 100. 

(RDYRT) = XH95 : The probability of being able to 

fill a demand at an arbitrary point in time multi- 
plied by 100. 

(ENOH) = XH94 : The average number of items on 

hand at an arbitrary point in time multiplied by 100. 

(ORD) = XH93 : The average number of orders 

placed in a unit period. 

(RI) = XH92 : The average number of repair 

inductions (batches of items repaired) per unit 
period. In the infinite server queueing system 
RI will be C. 

(C) = XH91 : The average number of items in 

the repair department at an arbitrary point in time 
multiplied by 100. 

(DIFF 1) = XH90 : The difference between ENOH and 

the theoretical value (V4) of the expected number of 
items on hand at an arbitrary point in time multi- 
plied by 100. 

(DIFF 2) = XH89 : The difference between B T and 

BARB multiplied by 100. This is calculated 1 to 
illustrate that the total time-weighted backorders 
incurred in a unit period is equivalent to the 
expected number of backorders at an arbitrary point 
in time. 
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NOTE: Measure of performance four (B„) would normally be 

normalized fF ) by dividing byten, i.e., 

XH97=1102= =1.102. 

the expected number of unit years of backorders and the 
average number of items on hand at an arbitrary point in 
time. For all runs the system was preloaded such that the 
inventory position was equal to the sum of the order quantity 
and reorder point. It was felt that this would result in 
negligible transient times, hence the entire run’s statistics 
were usable in all cases. The input values of all parameters, 
for what was called the parent system, are summarized in 
Table II. Deviations from these values of the parent 
system for any one experiment are pointed out in the . 
description of each experiment as well as in Table III. 



TABLE II. PARENT SYSTEM PARAMETER VALUES 



PARAMETER 



VALUE 



Order quantity (Q) 

Reorder point (R) 

Minimum batch size (L) 

Maximum batch size (K) 

Probability of repair (P) 

Mean demand during repair time (X) 
Type of distribution for repair time 
Mean repair time (X ) 



‘R- 



= 10 




= 0 



= 1.4 



=336 hours 



-constant 



Lower bound for uniform repair time (LBR) =0 
Upper bound for uniform repair time (UBR) =0 
Variance for repair time distribution «r R )=° 
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Type of distribution for lead time 




-constant 


Mean lead time (X^) 




= 504 


Lower bound for uniform lead time 


(LBL) 


= 0 


Upper bound for uniform lead time 


(UBR) 


= 0 


Variance for lead time distribution (6 "l) 


= 0 


Type of output desired 




-summary 



TABLE III. CHANGES TO PARENT SYSTEM PARAMETER VALUES 



A. CHANGES FOR EXPERIMENT ONE 
(Q)=l,2,3,...,10 

B. CHANGES FOR EXPERIMENT TWO 



(X) =1 . 4 , 2.75, 4, 5.5 



Type of distribution for repair time- -constant , 

uniform, normal, 
and empirical 

Type of distribution for lead time --constant, 

uniform, normal, 
and empirical 



LBR=168 hours 
UBR=504 hours 
6*R = 75 hours 



LBL=240 

UBL=768 

5“l =116 



C. CHANGES FOR EXPERIMENT THREE 

(Q)=2 00=1, 2, 3, 4, S, 6 

D. CHANGES FOR EXPERIMENT FOUR 

(Q)=2 (L)=l, 2, 3, 4 

E. CHANGES FOR EXPERIMENT FIVE 

(Q) = 2 (X)=4 

(X) =0 .40 , 0.65, 0.75, 0.80, 0.90, 1.15, 1.40 



00=1 

(K) =4 

(L) =l 
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B. DESCRIPTION OF EXPERIMENTS 

1 . First Experiment: The Whole Spectrum 

The first experiment considered the parent system 
but with values of mean demand during repair time of 1.4, 
2.75, 4, and 5.5. Also, order quantities of one through 
ten and variable reorder points were allowed. That is, for 
each combination of demand and order quantity above, the 
reorder point (R) was allowed to vary from an initial value 
of zero until the expected number of unit years of backorders 
reached an acceptable stationary level. This was done in 
order to see how the measures of performance, particularly 
the two mentioned above, varied with respect to changes in 
the order quantity and reorder point for the system described 
above . 

The amount of computer time required to run one period 
ranged from a low value of 5.89 seconds for (X)=1.4, with 
Q=10 and R=7 ; to a high value of 28.57 seconds for (X)=5.5, 
with Q=l, and R=0. The total computer time required for 
this experiment would obviously depend on the particular 
inventory system being investigated by the user. It is 
assumed that the user would know, determine, or estimate 
his system demand such that he could restrict his investiga- 
tion to one or two values. Therefore, the number of runs 
and total computer time required for his experiments would 
be significantly less than that required for these presented 
in this paper. Some typical results of the first experiment 
are shown in Figures 2 and 3. 
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2 . Second Experiment: Distribution Test 

For this experiment a representative system from 
experiment one was chosen for further study. The system 
chosen was thought to be a typical, low demand, high cost, 
repairable item inventory system with a reorder quantity of 
ten and a mean demand during repair time of 5.5. The par- 
ticular reorder quantity used is a somewhat subjective 
choice which the manager would make. Hopefully the manager 
could look at the results of experiment one and pick what 
would be the best quantity. A procedure to do this is 
explained in the analysis of experiment one. However, what 
will often happen is that one of several possible factors 
will force the manager to live with some predetermined order 
level. Some examples of these may be: there is a minimum 

number of items that the manager can let a contract for with 
a particular manufacturer, or an engineering requirement 
forces items to be procured in matched pairs or sets, or 
possibly a fixed unit of issue by the manufacturer. Regard- 
less of how it is decided upon, once the reorder quantity is 
fixed the sensitivity of the system can be investigated by 
changing various input factors such as the repair time 
distribution or lead time distribution, as was done in this 
experiment. Constant, uniform, normal, and empirical repair 
and lead time distributions were examined. For each combina- 
tion of distributions the reorder point (R) was varied until 
the total unit years of backorders was stationary. The 
purpose of this experiment was to show the sensitivity or 
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insensitivity of the modified parent system to various 
types of distributions for the repair and lead times each 
having the same mean and variance. Some typical results of 
the experiment are given in Figure 4. 

3 . Third experiment: Maximum Batch Size Test 

For this experiment the parent system was modified 
in order to treat^the repair department as a single server 
queueing system where the server is capable of handling as 
many as X items at one time. Using the procedure explained 
in the analysis of experiment one, the order quantity was 
fixed at a value of two for this experiment. Then the sen- 
sitivity of the reorder point (R) was investigated by varying 
R for each feasible value of X, the maximum batch size, until 
the expected number of unit years of backorders reached an 
acceptable stationary level. A simple way to determine the 
maximum value for X is to look at the results obtained in 
experiment one for the comparable infinite server queueing 
system. The run in experiment one where all the input 
parameter values were identical to this one with the excep- 
tion of the parameters associated with the repair discipline 
would be the comparable infinite server queueing system. One 
of the items included in the summary output for Storage 
Statistics, FASTR, summarizes the statistics for the case 
in which the repair department is treated as an infinite 
server queue. Two of the particular statistics included 
there are the average number of items in the repair depart- 
ment and the maximum number of items that were ever in the 
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repair department at any time. These particular values were 
used in this experiment as a guide for determining the 
maximum value for K to be considered. The results of this 
experiment are shown in Figure 5. 

4 . Fourth experiment: Minimum Batch Size Test 

In this experiment the parent system was again 
modified to treat the repair department as a single server 
queueing system. Again, as in the last experiment, the 
sensitivity of the reorder point to the type of repair 
policy used was investigated. This time, however, K 
rather than L was held constant while R was varied for each 
value of L, the minimum batch size, until the expected 
number of unit years of backorders reached an acceptable 
stationary level. The values of L were varied from one to 
K. The value of K used was four, which was the maximum 
feasible quantity for K as determined by experiment three. 
Since the repair department was treated as a single server 
queueing system, the repair of one batch of items must be 
completed before repairs can begin on the next batch of 
items. The results of this experiment are presented in 
Figure 6. 

5 . Fifth Experiment: Mean Repair Time Test 

The modified parent system of experiment four was 
used for this experiment but with L=1 and with various mean 
demands during repair time (X) . The mean demand during 
repair time was allowed to vary from 1.40 to 0.40 by 
changing the mean time to repair from 336 hours to 96 hours 
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while keeping the mean time between failures constant at 
242 hours. The sensitivity of the reorder point was then 
investigated by again varying R for each value of X until 
the expected number of unit years of backorders per unit 
period reached a stationary level. A physical interpretation 
here could be a change of equipment (perhaps manual repair to 
automation) or a change in the budget for the repair depart- 
ment. An investigation of the results of this experiment 
suggested a comparison of the trade off between decreasing 

the mean repair time (X ) and increasing the maximum batch 

R 

size in the repair department. The results of experiment 
five are presented in Figure 7. 
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IV. ANALYSIS AND CONCLUSIONS 



A. ANALYSIS OF THE FIRST EXPERIMENT 

In all experiments the reorder quantity and reorder 
level were controlled. In addition, for each experiment 
one variable such as the maximum or minimum repair batch 
size was controlled while the remaining variables were 
determined from the simulation. In this experiment, for 
each reorder quantity and mean demand during repair time 
considered, the reorder level which reduced the expected 
number of unit years of backorders per unit period to an 
acceptable level was determined. To determine this reorder 
point a one-way analysis of variance was used. 

Such a test requires additivity, linearity, normality, 
independence, and homogeneous variances. The first two 
assumptions were felt to be satisfied since the expected 
number of unit years of backorders per unit period appears 
to be a linear sum of the true mean effect, the true effect 
of the treatment (the value of the reorder point chosen) , 
and the true effect of the errors observed. The primary 
error in the observations is due to the stream of random 
numbers. It seems to be a good assumption that these errors 
are normally distributed with mean zero and variances which 
are equal for each treatment. 

This particular test was used because it is quite robust, 
that is, the test is reliable under rather strong modifications 
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to the above assumptions. It can also be shown that trans- 
formations of the original data are justified. This was 
particularly important in this case because it was necessary 
to make a logarithmic transformation to the observations 
in order to obtain a linear model. Figure 8 illustrates 
this type of transformation. 

The null hypothesis used in this analysis of variance 
was that the expected number of unit years of backorders 
per unit period for each of the reorder points was the same 
for a fixed reorder quantity and mean demand during repair 
time. For the parameters considered in this paper the 
analysis of variance was reported for that case which 
appeared to have the highest probability of not rejecting 
the null hypothesis. Nine observations with different 
streams of random numbers for each run were made and the 
null hypothesis was tested at the .01 significance level. 

The results of the analysis of variance yielded a test 
statistic of 194, well in excess of the F-value (2.48) 
required to reject the null hypothesis. (For calculations 
see page 79.) 

Since the primary goal of the first experiment was to 
determine the ’’best” order quantity, reorder point combina- 
tion, Scheffe’s [Ref. 5] multiple comparison test was then 
used to find this combination. This procedure tests for 
equality of means between contrasts in a systematic fashion 
until a significant difference between means is found. The 
analysis of variance indicated that a significant difference 
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was present. Scheffe f s test was then used to determine 
the minimum value of the difference in means which statistic 
ally differentiated among the various values of the expected 
number of unit years of backorder per unit period. 

In order to apply Scheffe ! s test a (Q,R) pair yielding a 
value of the measure of performance which was considered 
acceptable was chosen. The first contrast looked at the 
differences in the measure of performance for this pair 
(Q,R^) and the pair (Q,R^-1) . If this contrast failed to 
reject the null hypothesis, that the expected number of 
unit years of backorders per unit period were equal at the 
.01 significance level, then the next contrast would examine 
the pairs (Q,R^) and (Q,R^-2) . This procedure continued 
until the null hypothesis was rejected. If this happened 
when the pairs (Q,R^) and (Q,R -j) were compared, the 
minimum value of the differences which was statistically 
significant was known to be less than or equal to the 
difference between the values of the measure of performance 
for these two pairs. For future reference denote R^-j by R^ 
The minimum difference required for significance can be 
solved for exactly by applying this procedure in reverse, 
that is, by setting the value of an arbitrary contrast 
equal to the test statistic required to reject the null 
hypothesis and then transforming this contrast into the 
difference required in the measure of performance for 
significance. For the example systems considered in this 
paper this procedure required the use of three contrasts 
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before significance was found between (Q,R) pairs (10,11) 
and (10,8). This contrast resulted in a test statistic of 
6.99 which was greater than the F-value (2.48). Hence, the 
null hypothesis was rejected. (For calculations see pages 
81-83). Further, forcing the test statistic to be equal to 
2.48 and solving for the required minimum difference between 
values of the measure of performance for significance 
resulted in a value of 0.096 unit years of backorders per 
unit period. Statistically, one cannot distinguish a 
difference between 0.096 and 0.00 for this measure of 
performance. For this particular example, which is summar- 
ized in Tables IV and V, this means that any position to 
the right of the dashed line in Table IV could be an 
acceptable operating point for a system when X, the mean 
demand during repair time, was 5.5 and the other system 
characteristics were as noted in Tables II and III. 

The procedure then used to determine the ’’best" operating 
point for this particular example was to list the pairs 
(Q,Rq) and observe the measures of performance listed in 
Table I. Notice from Table V that when one tries to 
identify the pair that minimizes both the average number of 
items on hand at an arbitrary point in time and the average 
number of procurements per unit period the (Q,Rq) pairs 
(2,10) and (3,10) are selected as the "best” operating 
points for this example. A simple comparison of holding 
cost and procurement cost could enable the user to further 
distinguish between these two choices. 
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TABLE IV. TYPICAL RESULTS FROM EXPERIMENT ONE WITH B SIGNIFICANCE LINE 
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In applying these results to other systems it was ob- 
served that the value of 0.096 for the expected number of 
unit years of backorders per unit period could also be used 
to determine the (Q,Rq) candidates for the "best" operating 
point when considering all values of mean demand during 
repair time less than or equal to 5.5. It was further 
observed that among these (Q,R^) candidates, most measures 
of performance would differ by negligible amounts with the 
exception of the average number of items on hand at an 
arbitrary point in time and the average number of procure- 
ments made per unit period. 

B. ANALYSIS OF EXPERIMENTS TWO THROUGH FIVE 

1 . General Comments 

The analysis of variance and the multiple comparison 
test could have been u’sed to obtain candidates for optimal 
operating points with other relevant parameters, such as 
the minimum or maximum repair batch sizes or the distribution 
of repair time being used as the basis. Although this 
analysis was not repeated for all experiments, a close 
observation of the data seemed to indicate that the minimum 
significant difference value would be nearly the same 
regardless which parameter was chosen for the basis. Use 
of this observation was made in the analysis of experiments 
two through five. 

2 . Repair Time Distribution (Second Experiment) 

The results summarized in Tables VI and VII showed 
that the type of repair time distribution had a negligible 
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TABLE VI. TYPICAL RESULTS FROM EXPERIMENT TWO WITH B„ SIGNIFICANCE LINE DRAWN IN. 
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EMPIRICAL 1423 99.35 8 .03 .000 100 8.69 14 1292 4.96 -0.02 .073 



effect on the expected number of unit years of backorders per 
unit period provided the mean repair times were the same. 

This result was expected since Palm's theorem indicates that 
the distribution of the number of items in repair depends 
only on the repair time distribution through its mean value. 
When the probability of repair was nearly one it was also 
observed that the type of lead time distribution had a 
negligible effect on system performance. A probable 
explanation for this result is that the average number of 
orders placed per unit period was insignificant compared to 
the number of repairs made. This can be seen in Table VII. 

3 . Maximum Batch Size Test (Third Experiment) 

In this experiment the effects of change in the 
maximum repair batch size were investigated. The results 
summarized in Tables VIII and IX indicate a vast improvement 
in system performance for increasing the maximum batch size 
from one to two but negligible improvement thereafter. The 
initial improvement can be explained by the queueing result 
which says that in order to ensure stationarity the maximum 
repair batch size must be large enough to guarantee that the 
mean time between failures multiplied by K must be greater 
than the mean repair time. An important question which was 
considered was exactly how large K should be. In order to 
provide an answer further work would have to be done to 
determine an explicit rule for finding the "best” value of K. 
However, it is felt that this rule would depend on the 
utilization factor of the repair department. For this 
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particular example, increasing K to obtain a utilization 
factor of less than .417 resulted in negligible improvement 
in system performance. 

4 . Minimum Batch Size Test (Fourth Experiment) 

This experiment considered the effects of change in 
the minimum batch size for a fixed value of K. The results 
summarized in Tables X and XI indicated that, as before, 
changes in the minimum batch size requirement had little 
effect in overall system performance with the exception of 
the average number of repair inductions. If the cost of a 
repair induction is large certain economies could be achieved 
by requiring a minimum batch size. This could be achieved 
without degrading system performance. Further investigation 
would be needed to generalize these results in order to 
obtain a rule for determining the "best’' value for the 
minimum batch size. 

5 . Mean Repair Time Test (Fifth Experiment) 

For the case in which both K and L are one, the 
effects of change in the mean repair time were examined. 

The results summarized in Tables XII and XIII indicated that, 
as expected, the magnitude of the mean repair time is highly 
significant. In particular, system effectiveness was inversely 
proportional to the mean repair time. 
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TABLE X. TYPICAL RESULTS FROM EXPERIMENT FOUR WITH B_ SIGNIFICANCE LINE DRAWN IN. 
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339 99.12 3 0.03 0 97.56 4.44 17 213 1.81 - 0.60 .003 



TABLE XII. TYPICAL RESULTS FROM EXPERIMENT FIVE WITH B SIGNIFICANCE LINE DRAWN IN. 
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C. CONCLUDING COMMENTS ABOUT THE ANALYSIS 

It is interesting to note the tradeoffs between various 
parameters. For instance, it seemed that one could deter- 
mine a reduction in mean repair time which would give an 
equivalent increase in system performance to that provided 
by an increase in the maximum repair batch size. An illus- 
tration of such a tradeoff for the example considered in 
this paper can be seen by looking at Table XIII for the 
following entries: minimum and maximum batch size of one 

with mean repair time of 216 versus minimum batch size of 
one, maximum batch size of four, and mean repair time of 
336 hours. 

Decisions such $s determining the maximum batch size or 
mean repair time must be made by the decision maker based on 
the feasibility of change and the cost of such a change. 

This simulation will provide the manager with each of the 
measures of performance needed to make such decisions. 

In addition to the above remarks, another item should be 
discussed. Recall that one of the objectives of this study 
was to demonstrate the analytical results obtained by 
Richards [Ref. 7] . Notice that in each of the tables which 
summarized the measures of performance for the various 
experiments, two of the measures of performance listed are 
f, DIFF 1" and "DIFF 2.” "DIFF 1” reflects the difference 
between the simulated and theoretical values for the average 
number of items on hand at an arbitrary point in time. 

"DIFF 2" reflects the difference between the expected number 
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of unit years of backorders per unit period and the average 
number of items on backorder at an arbitrary point in time. 
In observing these differences in Table V it was noted that 
the maximum value for M DIFF 1" was 0.39 and for "DIFF 2" 
it was .009. This represents a maximum percentage differ- 
ence of 5.1% for f, DIFF 1.” It is felt that since these 
differences are so small, this simulation does in fact 
demonstrate the validity of the analytical results obtained 
by Richards . 



V. PROGRAM GUIDE 



This section is written so that, along with Table I 
(Supply Performance Measures) and the table of symbols and 
abbreviations, it can be a self contained unit for the user. 
It contains all the information needed to successfully set 
up and run the model. 

A. GENERAL COMMENTS 

Most people will agree that any inventory system should 
be designed to give the maximum amount of service for the 
minimum amount of cost possible. Hence most inventory 
studies evaluate the various policies being considered by 
looking at the costs incurred over a unit period of time. 
Then policy A is chosen over policy B if the total expected 
cost per unit time using policy A is less than when using 
policy B. However, quite often in this approach one cannot 
specify all the cost elements. An alternative approach is 
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to determine the average values of such measures of perform- 
ance as the number of items on hand at an arbitrary point 
in time, the number of backorders per unit period, and others 
as listed in Table I. The user can then choose from these 
measures of performance those he wishes to use to evaluate 
the policies under his scrutiny. This method provides the 
user with an easy way of deciding which values for system 
parameters such as the reorder quantity and reorder point 
yield satisfactory performance. 

Closed form solutions can be obtained for the various 
measures of performance but often the solutions are not very 
tractable for further calculations. This is where the 
simulation is of value. It systematically and quickly 
calculates the magnitude of all the measures of performance 
for various input values of all parameters involved and 
presents the user with a fast and easy means of investigating 
the sensitivity of his inventory system to any changes in 
these parameters. To try and do this procedure by hand or 
apply the analytical results to a sensitivity analysis of 
this size would involve a large number of man-hours of work 
to produce effective results. 

B. EVENTS 

This model is an event-step computer simulation written 
in the GPSS/360 language, i.e., all actions that are to 
occur in the simulation are dynamically generated by the 
computer program as a result of previous simulation actions 
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and are listed chronologically in order of descending 
priority in a Current Events Chain, Future Events Chain, 
or Interrupt Chain. For purposes of programming, each of 
the major actions included in the simulation can be con- 
sidered as an entity in itself and assumes the form of a 
computer program subroutine. Such major actions have been 
called events for convenience . The eight major actions 
(events) included in the simulation are: 

(1) Item Failure 

(2) Item Inspection 

(3) Scrap 

(4) Repair 

(5) Issue 

(6) Backorder 

(7) Order 

(8) Receipt of Material 

Each of the computer program subroutines representing 
these events uses as input parameters the following informa 
tion: 

(1) Time the event is to occur 

(2) Identification of event 

(3) Number of items involved 

(4) Priority level (0,1,...) 

(5) Current values of system characteristics such as 
inventory position (IP) , number of items ready for 
issue (D) , and others. 

The dynamic process of simulating one inventory policy 
from start to finish forms the executive routine for the 
computer simulation. This executive routine also includes 
three program subroutines referred to as Random Observation 
Calculation of Measures of Performance, and Incrementation 
of System Parameters. General flow charts describing the 
logic included in each event of the simulation plus the 
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interrelationship of events are included as Figures 9 through 
12. Detailed flow charts for the entire model are provided 
in Appendix B. Random observation is the subroutine that 
generates random points in time at which certain parameter 
values are recorded in order to obtain certain useful sta- 
tistics about the system. Calculation of Measures of 
Performance is the subroutine that, at the completion of 
any unit period (a time period which represents 10 years 
of actual system operations) , interrogates the information 
contained in all system parameters and transfers control of 
the computer program to the calculation phase of the measures 
of performance listed in Table I. Incrementation of system 
parameters is the subroutine which increments the current 
value of R, the reorder point, by one. It also clears the 
results of the last unit period from the program and resets 
all the random number generators to the same values used in 
the previous run to ensure meaningful comparisons between 
runs . 

C. INPUT 

All inputs to the computer program are contained on 
three punched cards called Initial Cards and four Function 
Follower Cards. The format for this information and the 
units to be used are detailed in Tables XIV through XVI. 

With one set of input cards a user can simulate the inter- 
actions involved for as long as 140 years or 14 unit periods 
of ten years each. Each unit period will increment the 
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Figure 9. Relationships of events. 
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BEGIN (I ) ETEM FAILURE 



(^INITIALIZER) 
f ^GENERATE 



M = M + 



2 ) INSPECTION 



(2) ITEM INSPECTION 




TYPE 2 



(3) SCRAP (TYPE 2) 




BACKORDER 




/-ADVANCE^ 
( CLOCK BY ^ 
R EPAIR TIME 



LEAVE 
REPAIR DEPT 
C= C- 1 



BACK- 

ORDER 



ISSUE 



RECEIPT OF 
MATERIAL 



Figure 10. Major Events (I ) — ( 4) Logic Diagrams. 
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Figure II. Major Events 



7 ) (7) ORDER 




ORD=ORD+l 
A= A +Q 
IP = IP +Q 



f ENTER N 
B ORDER DEPJ 7 
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BY LEADTIME 



LEAVE 
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/TEST' 



RECEIPT OF 
MATERIAL 

(8) RECEIPT OF 
MATERIAL 



<^ce fI p>>no_0 



YES 



D = D +1 



ISSUE 




(5) — (8) Logic Diagrams. 
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(II) Incrementation of system parameters and control of program 
running time. 




This clears all values from program. 
It also resets the random number 
generator so the same number of 
XAs are created for each succes- 
sive unit period. 



Figure 12. 

Administrative events (9) -(II) logic 
diagrams. 
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value of R by one starting with the initial value chosen by 
the user (generally zero) . In order for the user to vary 
parameter values other than the reorder level, a new set of 
initial cards must be prepared by the user for the program. 

All input values must be integers and in units of hours if 
it is a time input. For example in the three columns allotted 
for Q on card one the user could enter anything from 001 to 
999 for one item up to 999 items or in the columns allotted 
for the mean time between demands (failures) , on card 
two the user could enter anything from 999 to 001 to repre- 
sent failures rates of .7308 per month (one every 999 hours) 
to 730 per month (one every hour) . 



D. FORMAT SECTION 

This section contains the tables which identify and 
explain the various system parameters and characteristics 
needed as user input in order to use the simulation. It 
also includes the feasible range for each entity. This 
input must be typed by the user on three Initial Cards and 
four Function Follower Cards. 

1 . Description Of Entities Found On The Initial Cards 
This input is summarized in Tables XIV through XVI. 

2 . Description Of Function Follower Cards One and Two 
These cards contain the function values, Y^, and 

argument values, X^ , for the empirical distribution of 
repair time. This is a discrete numerical valued function 
with up to ten different points. Figure 12 shows that this 
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TABLE XIV. FIRST INITIAL CARD 



Q Reorder quantity. Number of items reordered at 

one time. 

R Reorder Point. If Inventory Position is < R an 

order is placed. 

L Minimum Batch Size . Under the batch repair policy, 

repair won’t start unless at least L units are 
currently waiting to be repaired. Ensure L=0 if 
the infinite server queueing policy for repair is 
desired. 

K Maximum Batch Size. Under the batch repair policy, 

this is the maximum number of items that can be 
repaired at one time. 

P Probability of Repair. The probability that a 

failed item can be repaired. If P=0000, (P=0.0) 
the model degenerates to the consumables case. 

If P=1000, (P=1.0) the model becomes one without 
degredation of supply which implies that no orders 
are placed and no outside procurement is needed. 

.. Mean time between failures. The expected number of 

'u hours between demand. 

FLAG/R Repair Time Distribution Flag. One of four numbers. 

3 Constant repair time. 

4 Uniform distribution for repair time. 

5 Normal distribution for repair time. 

6 Empirical distribution for repair time. 
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TABLE XV. SECOND INITIAL CARD 




Mean Time For Repair. The expected time required 
for repair of an item or batch of items. 


LBR 


Minimum Repair Time. The lower bound on repair 
time for the uniform distribution of repair time 
(used when FLAG/R=4) . 


UBR 


Maximum Repair Time. The upper bound on repair 
time for the uniform distribution of repair time, 
(used when FLAG/R=4) . 




Standard Deviation For Repair Time. The standard 
deviation for the normal distribution of repair 
time, (used when FLAG/R=5) . 


FLAG/R 


Lead Time Distribution Flag. One of four numbers. 

7 Constant lead time. 

8 Uniform Distribution for lead time. 

9 Normal Distribution for lead time. 

10 Empirical Distribution for lead time. 




Mean Time For Procurement. The expected value for 
procurement lead time. 

TABLE XVI. THIRD INITIAL CARD 


LBL 


Minimum Lead Time. The lower bound on lead time 
for the uniform distribution of lead time. 


UBL 


Maximum Lead Time. The upper bound on lead time 
for the uniform distribution of lead time. 


*L 


Standard Deviation For Lead Time. The standard 
deviation for the normal distribution of lead time. 


XH17 


Type Of Output Desired Flag. One of two numbers. 
Zero for summary output, or one for standard output. 
(Normally a zero should be used to conserve printing 
time . ) 
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function has the same function value, FN.=Y., for all 

3 i 

argument values in the interval between X.-l and X.. No 

J J 

interpolation is performed, and the vlaue at the right- 

hand end of the interval is used. For argument values 

less than the magnitude of X^ , FN^=Y^. For argument values 

greater than the last value of X., FN =Y . 

& j 9 n n 

The "function follower cards contain the individual 
points (successive pairs of X^. , Y values) of the function. 
For example, consider the function pictured in Figure 12 as 
the repair time distribution, the function follower cards 
would then be: 0.06,439/0.07,439/0.08,439/0.58,464/0.66,486/ 

0.68,558/0.69,558/0.71,558/0.96,582/1.0,606. 



Notice that in this example there are really only 
six different function values FN^ , so argument values 
.06, .07, .68, and .69 are dummy values and are put in only 
to meet the requirement of having ten pairs of X^. ,Y values. 

3 . Description Of Function Follower Cards Three and Four 
These cards are used exactly the same way as Function 
Follower Cards One and Two except they contain the function 
and argument values for the empirical distribution of lead 
time. For example, if the function pictured in Figure 13 
were the lead time distribution, FFC’s three and four would 
be: 0.2,360/0.4,384/0.5,432/0.6,480/0.8,768/0.9,840/ 

0.95,864/0.99,888/0.995,912/0.999,984. 



Notice that once FLAG/R and FLAG/L have been 
selected the program ignores those cards which involve the 
parameters of the other distributions. Therefore only those 
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V, = Constant 



FN 




Discrete numerical valued function with 10 arguments and duplicate 
function values for sample empirical repair time. 




Discrete numerical valued function with 10 arguments, no duplicate 
function values but with argument X n < 1.0. \ function values for 
.TI5<X<I.O FN =984. (sample empirical lead time.) 
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values or cards which affect the parameter values of the 
particular type of distribution specified need be worried 
about. For example, if FLAG/R=5 and FLAG/L=10, the user 
must ensure and have the values desired and FCC's 

three and four have the values desired. The values of LBR, 
UBR, LBL, UBL , 5^, and FFC ' s one and two can be ignored. 

4. Examples of Initial Cards 

Tables XVII through XIX give examples of the three 
Initial Cards. 

TABLE XVII. EXAMPLE OF INITIAL CARD ONE 
Columns Name Range Units 

8-14 INITIAL 



19-22 


xhi; 






23-24 


Q 


1-99 


ITEMS 


25-29 


/XH2 , 






30-31 


R 


0-99 


ITEMS 


32-36 


/XH3 , 






37-38 


L 


0-K 


ITEMS 


39-43 


/XH4, 






44-45 


K 


0-99 


ITEMS 


46-50 


/XHS , 






51-54 


P 


0-1000 


.001 


55-59 


/XH6 , 






60-62 


-^0 


0-999 


HOURS 


63-67 


/XH7, 






68-69 


FLAG/R 


3,4,5, or 6 


NONE 



All input information is right- j ustified in the 
space allotted. Example: 

COLUMN NUMBER - 28 39 55 69 

tt ♦ t 

INITIAL XH1 , 10/XH2 .01/XH3 , 00/XH4 , 01/XH5 , 0900/XH6 ,061/XH7,03 
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TABLE XVIII. EXAMPLE OF INITIAL CARD 


TWO 


Columns 


Name 


Range 


Units 


8-14 


INITIAL 






19-22 

23-25 


XH8 , 


0-999 


HOURS 


26-30 

31-33 


/XH9, 

LBR 


0-998 


HOURS 


34-39 

40-43 


/XH10, 

UBR 


0-998 


HOURS 


44-49 

50-52 


/XH11, 

«k 


0-580 


HOURS 


53-58 

59-60 


/XH12 , 
FLAG/L 


7,8,9, or 10 


NONE 


61-66 

67-69 


/XH13 , 


0-999 


HOURS 


All input information 
space allotted. Example: 


is right -justified in the 


COLUMN 

INITIAL 


NUMBER - 34 

t 

XH8 , 336/SH9 , 168/XH10 , 


53 

t 

0504/XH11 , 095/XH12 , 07/XH13 




TABLE XIX. EXAMPLE OF 


INITIAL CARD THREE 


Columns 


Name 


Range 


Units 


8-14 


INITIAL 






19-23 

24-26 


XH14, 

LBL 


0-998 


HOURS 


27-32 

33-36 


/XH15, 

UBL 


0-1998 


HOURS 


37-42 

43-45 


/XH16, 

®L 


0-580 


HOURS 


46-51 

52 


/XH17 , 


0 , or 1 


NONE 


space 


All input information 
allotted. Example: 


is right- justified in the 


COLUMN NUMBER - 2 6 37 51 

t * t 

INITIAL XH14 , 240/XH15 ,0768/XH16 , 116/XH17 ,0 





65 



E . OUTPUT 



The computer program has two different sets of output, 
summary (special) and standard. The type of output received 
depends on the value of XH17. If XH17=0 the summary output 
is received, but if XH17=1 the standard output is received. 
Computer running and printing time is greatly increased if 
the standard output is used. Therefore, it is recommended 
that XH17 be set to zero. If there is a major error in the 
program, the standard output will be printed out automatic- 
ally. Following is an explanation of the two sets of output. 

1 . Summary Output 

An example of the summary output is illustrated in 
Appendix D. Notice that there are three major divisions: 

(1) savevalues, (2) storage statistics and (3) queue 
statistics . 

a. Savevalues 

There are four subdivisions in this part. The 
first is a list of nonzero halfword savevalues between 
five and 100 for system parameters and statistics. This 
list is included for completeness only and not because of 
their nature or importance. The second subdivision, Clock 
Time, merely shows the number of hours of operation 
simulated. Normally this will be 87600 hours or ten years 
of operation, which has been defined to be a unit period. 

The third subdivision, Critical System Parameters, lists 
the values of Q, K, and L which the user has chosen and 
input on the initial cards. It also prints out the value 
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of R for this particular time period. The last subdivision, 
Measures of Performance, includes most of the important 
statistics to be considered when judging or comparing 
alternative inventory systems. 

b. Storage Statistics 

This division provides the important statistics 
for the order department (ORDER) , the infinite server 
queueing repair process (FASTR) , and the single server 
queueing repair process (SLOWR) . The important items and 
explanation of each storage follow. (1) Storage -- the name 
used to represent the physical entity. (2) Capacity -- the 
maximum number of transactions which can use the storage 
at any one time. (3) Average Contents -- the average 
number of transactions that were in the storage during the 
run. (4) Entries -- the number of transactions which 
entered the storage. For SLOWR this is the number of 
batches of items that entered the repair process rather 
than the number that were repaired as in the FASTR storage. 
Likewise, it is the number of orders placed rather than the 
number of items ordered for the ORDER storage. (5) Average 
Time/Tran -- the average number of hours that an order was 
outstanding, or the average number of hours that it took 
for an item to be repaired, or the average number of hours 
it took for a batch of items to be repaired. (6) Maximum 
Contents -- the maximum number of transactions that used 
the storage at any time during this run. 
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c. Queue Statistics 

This division reflects the important statistics 
for the infinite server queue repair department (FSTRQ) and 
the single server queue repair department (SLOWQ) . The 
important items and explanation of each follow. (1) Queue - 
the name used to represent the physical entity. (2) The 
explanation for maximum contents, average contents, and 
average time/tran are the same as for storages. (3) Total 
Entries -- the number of transactions which entered the 
queue. (4) Zero Entries -- the number of transactions 
which entered the queue and departed with no time delay. 

Note that for the FSTRQ it should always be the same as the 
total entries. (5) Average Time/Trans -- the average number 
of hours that the delayed transactions, only, spent waiting 
to be repaired. 

2 . Standard Output 

An example of the standard output is illustrated in 
Appendix D. Notice that the standard output includes all of 
the items listed in the summary output as well as the follow 
ing items. (1) table statistics, (2) group counts, and 
(3) block counts. 

The table statistics provide the statistics of the 
following tables: BTCHZ , BTIME , ENOH, NBARB , NUMC , and 

RDYRT. The important items and explanations of each follow. 
(1) Table -- the name used to collect the particular 
statistic. (2) Entries in Table -- the number of times 
that an entry was made in the table. (3) Mean Argument -- 
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the average value of the entries made in the table. 

(4) Standard Deviation -- the standard deviation of the 
entries made in the table. (5) Upper Limit -- the maximum 
value of each cell of the table. For example, all entries 
between 51 and 100 fall in cell two of BTIME table. 

(6) Observed Frequency -- the number of entries which fall 
into this cell. (7) Percentage of Total and Cumulative 
Percentage are self-explanatory. (8) Multiple of Mean -- 
the upper limit of the cell expressed as a multiple of the 
mean of the table. The additional pieces of information 
given in the tables as well as the queue and storage 
statistics which are not explained in detail here are those 
which, under normal circumstances, are of no concern to the 
user. However, there may be occasions when a user wishes 
to make changes internal to the program such as adding the 
capability of handling compound Poisson demands. In such 
a case the standard output would be a necessity, especially 
the block counts, in helping the user debug the program. 

The reason this output was included is to show the user 
that this additional information is available from the 
computer, and to show the user the type of output he will 
receive by default should an error be made while using the 
model. If the user desires to receive the standard output, 
XH17 should be set to one on initial card three. If one 
desires further information relating to computer output, 
he should see the IBM User's Manual [Ref. 4]. 
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VI. USES AND APPLICATIONS OF THE MODEL 



This model was built for three major reasons. First 
and foremost was to give a convenient and useful tool to 
present day inventory managers faced with the decisions of 
how much to buy, when to buy, and how to conduct repairs 
for the items under their cognizance. As was pointed out 
earlier, there are presently two major approaches in vogue 
to help the manager solve his problem; cost analysis and the 
use of simple procurement policies with specified measures 
of supply performance. Under the cost analysis approach 
the user runs into the difficult problems of quantifying 
costs. For instance I, the inventory holding cost, consists 
of many things, some of which all managers agree on, such 
as storage and handling, and some which they don’t agree 
on, such as loss, breakage, and opportunity cost. Similar 
disagreements concern the fixed order cost as well. For 
example, the Navy uses a sliding scale from $25 to $75, 
while the Army uses a different sliding scale from $50 to 
$400. Further problems arise when one tries to assign a 
cost to a critical item being out of stock when needed. 
Primarily because of these difficulties in quantifying 
costs, the second approach was used in this paper. It 
was felt that statements such as ”a fill rate of 951 must 
be maintained" would at least mean the same thing to all 
managers as well as being something that all managers could 
agree on. 
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The second reason this model was built was to provide 
an efficient way of obtaining insight into the problems 
faced by the inventory manager. Simulation is a straight- 
forward and logical way of doing this. In addition simula- 
tion avoids some of the pitfalls of the more classical 
means of analysis. For example, the model can use any 
time-step desired and simulate as far into the future as 
desired. The model easily handled the different types of 
repair processes which analytically prove to be quite 
cumbersome. When these reasons are coupled with the 
combined effects of uncertainty, dynamic interaction of 
system events, and interdependence among system parameters, 
it is easily seen why this problem was particularly amenable 
to simulation. Richards [Ref. 7] showed that analytical 
expressions could be obtained for the individual measures 
of performance used in this paper. However, to actually 
generate numbers from these expressions for a sensitivity 
analysis of this size would be a monumental task. Therefore, 
simulation again seemed to be called for. Also, the advan- 
tage of time compression, both in being able to simulate 
into the future in a matter of seconds and being able to 
obtain numerical values for the various measures of per- 
formance, made simulation appealing. 

A major criticism of simulation is that the simulated 
queueing model yields only an estimate of parameters, 
such as a waiting line’s average length or associated 
probability of delay, and as such is subject to statistical 
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error. A second criticism is that if the system is so 
complicated and complex that the other methods fail and a 
simulation is required, it will necessarily be a very 
difficult job to build the model and analyze its results. 

In this particular case the first disadvantage acted 
to a much lesser degree since analytical expressions for 
all the measures of performance used in this model do exist. 
This fact allowed less replication to obtain trustworthy 
answers in the search for the best choices of the policy 
parameters. A moderate number of runs in a reasonable 
period of time accompanied with the statistical tests of 
section four also lent credibility to the answers and 
conclusions obtained. The second disadvantage noted went 
for naught since the model exists and works. 

Finally, the third reason that this model was constructed 
was to demonstrate the validity of certain analytical 
expressions for some of the measures of performance used 
in this paper. As was pointed out in section four, the 
model succeeded in doing this. From the comments above 
and the results presented in this paper, it is easily seen 
that the model can be used effectively to describe current 
inventory systems, to explore hypothetical inventory 
systems, and to design or improve given inventory systems, 
as well as be used as a teaching aid. The simulation 
model can be used to help answer such a current question 
as how to best conduct repairs on the engines for the new 
class of destroyers DD963. 
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In applying this model to any of the uses mentioned 
above, one is always faced with the problem of deciding 
which measure of performance to use as the criterion for 
making decisions. In this paper the expected number of 
unit years of backorders per unit period was used. One 
reason it was used over the average number of items on 
backorder per unit period was that it considers both the 
number of backorders and the amount of time the backorders 
were outstanding, while the other merely considers the 
number of backorders, totally ignoring the amount of time 
the backorders were outstanding. It was felt that both 
the fill rate and the ready rate were even worse than the 
average number of backorders since they don't even take 
into consideration the number of units backordered, but 
merely note whether or not any backorders exist. Perhaps 
the most significant reason the expected number of unit 
years of backorders per unit period was used as the criterion 
was that once it was considered acceptable, the remaining 
measures of performance were also acceptable. 

There were occasions when another measure of performance, 
the average number of items on hand at an arbitrary point 
in time, was used to help determine an optimal operating 
point for a given system. Such was the case in experiment 
one when there appeared to be a number of acceptable 
operating points, (Q,R) pairs, for the system under investi- 
gation. This measure of performance was used in these 
crises to differentiate between these pairs. 
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As mentioned earlier, this model can also be used as a 
teaching aid as well as demonstrating analytical results 
obtained by others. Table XX (page 75) is a good example 
of this. If one sets the probability of repair equal to 
one, the model degenerates to the all repairables case 
without degradation. From the table, one can see how the 
various measures of performance -.change as the mean demand 
during repair time varies. Further, if Q is required to 
be one, the results of this model can be compared with the 
results of Feeney and Sherbrooke [Ref. 6]. Notice that in 
their results the fill rate was always zero when the reorder 
point was zero. This result was a consequence of their 
initial conditions, whereas the results obtained by this 
model do not include that restriction. If we again restrict 
the model by considering the all consumables case (i.e. 
letting the probability of repair be equal to zero) the 
model becomes that described by Hadley and Whitten [Ref. 3]. 
Using a sample problem, one again has an example of how the 
model can be used to verify known analytical results. This 
example and comparison to the cost approach is summarized 
in Table XXI. 

Of course there are endless experiments one could carry 
out with the use of this model. For instance one could use 
this model to help decide which items should be designated 
repairable. This could be done by comparing the measures 
of performance obtained by considering an item first as a 
consumable and then as a repairable. It would also be easy 
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TABLE XX. COMPARISON BETWEEN FEENEY § SHERBROOKE* S [REF. 6] RESULTS WITH THIS 
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to extend this model to make use of the variance reduction 
techniques explained in reference nine or to add the 
capability of handling compound Poisson demands. 

VII. SUMMARY AND CONCLUSION 

The analysis of the repairable inventory problem pre- 
sented in this paper was based upon a simulation model 
thought to be sufficiently realistic and general enough to 
represent many different inventory systems. 

The results have shown that the simulation can be a 
helpful tool for the managers of inventory systems to use 
to make wise decisions governing their policies concerning 
when and how much to reorder. It was demonstrated in this 
paper that the mean time to repair an item, but not the 
repair time distribution itself, was extremely critical to 
system performance. It was demonstrated, at least for 
systems with a mean repair time less than or equal to 1.4, 
that a constraint placed on the minimum repair batch size 
had a negligible effect on system performance. It was 
further demonstrated for this same system that reducing 
the mean repair time by 401 was more profitable than doubling 
the maximum repair batch size. It also demonstrated that 
increasing maximum repair batch size beyond those values 
which gave a service rate of .417 or less had negligible 
effect in terms of increasing system performance. This 
paper also emphasized the power of using the expected 
number of unit years of backorders per unit period as the 
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criterion for judging system performance. This measure, 
when reduced to an acceptable level, was accompanied by 
good values of all other measures of performance as well. 
This is not true of the other measures. For example, if 
the fill rate had been used, Table XIII shows that an 
acceptable fill rate doesn't necessarily imply that all 
the other measures of performance are also acceptable. 

It is believed that this paper supplements the tools 
which are needed to make intelligent decisions concerning 
the repairable item inventory system. Of course it can 
also be used for the consumables case by merely setting the 
probability of repair equal to zero. 

A couple of further remarks seem to be in order. First, 
it would seem worthwhile to determine if current repair 
departments could be utilized more fully. Perhaps this 
could be achieved by spending a greater fraction of budget 
dollars to repair available carcasses rather than continuing 
to procure new material. Secondly, this study seems to 
indicate that study of a stock point's physical plant for 
repairables could be profitable. For example, the fourth 
experiment indicated that the minimum repair batch size 
required before repair begins on a batch of items could 
fluctuate to some degree before system performance is 
degredated. 

Much has been learned about the interrelationships and 
effects of many of the parameters of repairable and non- 
repayable inventory systems. Perhaps this has been the 
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most significant result of the project. It would seem 
impossible for anyone who would work with this simulation 
not to come away with a better understanding of the hidden 
mechanisms and inner workings of an inventory system. 

In closing, a simulation model for continuous review 
inventory systems was presented. The model used, demonstrat- 
ed the sensitivity of Q, R, K, and L to numerous operational 
constraints and system parameters. A class of simple 
policies with various measures of performance was used 
rather than a traditional cost postulated model to determine 
acceptable parameters for a given inventory system. It is 
felt that the model can be successfully implemented in a 
very large number of inventory problems. 
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APPENDIX A. STATISTICAL TESTS 



Analysis of variance. This test was used to find the 
significant difference in observed values of the expected 
number of unit years of backorders per unit period. 





Sum of 
Squares 


Degrees of 
Freedom 


Mean 

Square 


F Ratio 


Between groups 


436.3584 


11 


39.6689 


194.0706 


Within groups 


18.8052 


92 


0.2044 




Total 


455.1636 


103 







Table value from' F-Distribution for 92 and 11 degrees of 
freedom at .01 significance level = 2.48. 

194.0706^ 2.48, therefore the null hypothesis was 

rej ected. 

Null hypothesis: The observed values for the expected 

number of unit years of backorders per unit period are the 
same for each treatment (reorder point value) . 
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<><> 



K. A. Scheffe’s Multiple Comparison Test. This test 
was used to find the magnitude required for the expected 
number of unit years of backorders per unit period to be 
significantly different from zero. 



If 



(i c ,i 



(N-1)(V'(C.) ) 



> 



(l-°0 ( n 1 ) (n 2 )‘ 



then the hypothesis that the contrast C_. differs from 
zero is rejected. 



C. = 
1 



t 

5 



C. . T. 
il i 



( V 

(T i ) 



T. 

l 



c ij 



N-l 

(n^Cnp 



is the estimated variance of contrast . 

t h 2 

is the estimated variance of the i tn treatment=S . 

is the sum of the n. observations for the i^ 
i 

treatment . 

is the mean square of the experimental error, 
is the number of treatments included in contrast C^. . 
are the coefficients for the j*'* 1 comparison with 
the i^* 1 treatment. 

is the significance level chosen=.01 for this 
particular test. 

is the total number of treatments less one. 
are the degrees of freedom. 
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Let be the contrast between (QR) pairs (10,11) and 
(10,10) respectively. Then: 

C 1= C 10 i T 10 + C 11>1 T 11 =(3)(-25.3176)+(-4)(-22.1442)=12.62 

0(0,)= X C 2 (T.)=(3 2 ) C8) (.2044)+(-4) 2 (6) (. 2044) =34 . 339 
1 i=10 il 1 

(|c n |) 2 = (12. 62 ) 2 = .422 

(N-ijttf (11)(34.339) 

.422 2.48 therefore, H q is not rejected. 

Let C 2 be the contrast between (QR) pairs (10,11) and 
(10,9) respectively. Then: 

C 2~ C g 2 T 9 +C 11,2 T 11 =(4:) (22. 404) + (-6) (- 22 . 1442) =43 . 24 

V^ C 2-* = (4 2 ) (9) (.2044) + (-6) 2 (6) (.2044) =73. 584 

( lc,l ) 2 = (43 . 24) 2 = 2.309 2.309 2.48 

(N-l) ((7(C 2 )) (11) (73.584) Therefore , H Q 

is not 
re j ected . 

Let be the contrast between (QR) pairs (10,11) and 
(10,8) respectively. Then: 



C 3 =C 8 3 T 8 +C 11,3 T 11 = ^ 4 ^ C' 14 • 35 3) + (-6) (-22 . 1442) =75 . 453 
V (Cj) = (4) 2 (9) (.2044) + (-6) 2 (6) (.2044) =73. 584 



3 = (75 .453) 2 = 6.994 6.994 > 2.48 

(N-l) ((? (C,)) (11) (73.584) Therefore, 

reject H q , the 
null hypothesis 
that the means 
are equal . 
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Let be the contrast between (QR) pairs (10,11) and 
(10, m) respectively. Then: 

= 2.48 = C| C m | ) 2 = (2.48) (11) (73.584) = 

2007.372 



C ! C m| ) 2 

(N-l)(A(C m J) 



C m = 44.8 (minimum contrast value 

required for statistical 
significance . ) 

Then: 44.8=C x>m T x+ C 11>mTll =(4)T x+ (-6) (22.404) 



T x = 22.016 
log e B t = -2.450 

B = .086 the significant difference of the expected number 
of unit years of backorders from .01 

B t = . 01+ . 086=0 . 096 the minimum significant difference 
from zero 
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APPENDIX B 



Detailed Flow Charts 
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THIS IS A TYPE2 ITEM (NON REPAIRABLE) SEE IF WE CAN FILL DEMAND 
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STATISTICS FOR REPAIR DEPT (FINITE SYSTEM) 
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